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ABSTRACT 



The quasilinear mean-field theory for driven MHD turbulence leads to the result that the observed cross helicity (it ■ b) 
may directly yield the magnetic eddy diffusivity r/T of the quiet Sun. In order to model the cross helicity at the solar 
surface, magnetoconvection under the presence of a vertical large-scale magnetic field is simulated with the nonlinear 
MHD code NIRVANA. The very robust result of the calculations is that (u z b z ) ~ 2{u • b) independent of the applied 
magnetic field amplitude. The correlation coefficient for the cross helicity is about 10%. Of similar robustness is the 
finding that the rms value of the magnetic perturbations exceeds the mean-field amplitude (only) by a factor of five. 
The characteristic helicity speed u v as the ratio of the eddy diffusivity and the density scale height for an isothermal 
sound velocity of 6.6 km/s proves to be 1.1 km/s. This value well coincides with empirical results obtained from the data 
of the HINODE satellite and the Swedish 1-m Solar Telescope (SST) providing the cross helicity component (u z b z ). 
Both simulations and observations thus lead to a numerical value of r]T — 1.1 X 10 12 cm 2 /s as characteristic for the 
surface of the quiet Sun. 
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1. Introduction 

It is not easy to measure the turbulent magnetic diffusiv- 
ity ffr at the solar surface. This quantity determines the 
decay of magnetic magnetic structures with scales larger 
than those of the turbulence. Theoretically, the decay of 
the magnetic structures should depend on the relation of 
the magnetic field amplitude to the so-called equipartition 
value B oq = J fj,Qp(u 2 ) defined by the turbulence. This phe- 
nomenon is known as the effect of 77-qucnching, i.e. the sup- 
pression of the eddy diffusivity by the magnetic field. 

The simplest realization of the ?y-quenching at the so- 
lar surface can be given with two numbers. The decay of 
active regions after Schrijver & Martin (1990) can be un- 
derstood with an eddy diffusivity of 10 12 cm 2 /s while the 
decay of sunspots which their much stronger fields leads 
to 10 11 cm 2 /s (Stix 1989). These values are smaller than 
the value of 3 • 10 12 cm 2 /s which results from the widely 



culated atmosphere models then fluctuation measurements 
can be used to calculate the numerical value of the eddy 
diffusivity for weak fields. This the more as (u ■ b) ~ (u z b z ) 
if the large-scale magnetic field has only a vertical compo- 
nent and the only vertical gradient is due to the density 
stratification. The correlation of the vertical components of 
flow and field can empirically be obtained by both Doppler 
measurements and spectropolarimetry. 

To estimate the value of the cross helicity we assume a 
density scale height of 100 km and write the result in the 
form 



("■*») 

B z 



H 7 



km/s, 



(1) 



used formula t]t 



c r ,u r 



ms^corr 



with the tuning parameter 



c v ~ 0.3, the correlation length £ corr and parameter val- 
ues taken close to the surface. Up to now there was no 
possibility to observe the turbulent diffusivity on the solar 
surface for the quiet Sun where the magnetic quenching of 
this quantity by large-scale magnetic fields is negligible. 

Riidiger ct al. (2011) have shown that the combination 
of a vertical field with a driven turbulence in a density strat- 
ified medium leads to an anticorrclation of the cross helicity 
and the vertical large-scale field, i.e. (u ■ b) = —rfiB z / H p 
with Hp as the scale height of the density. If both the cross 
helicity and the large-scale vertical field are known then 
the ratio of the eddy diffusivity and the density scale can 
be computed. If also the density scale is known from cal- 



wherc H7 — H p /100 km and tit — 10 7712 cm 2 /s. With ob- 
servations of the LHS of (TIJ of about 1 km/s one would find 
t?t of order 10 12 cm 2 /s. In the present paper numerical sim- 
ulations of stratified magnetoconvection and observational 
results are discussed and the theory will also be extended 
to the inclusion of a vertical stratification of the turbulence 
intensity. Both the simulations as well as the observations 
lead to very similar results for the desired magnetic eddy 
diffusivity for the quiet Sun exceeding the value given by 
(HJ) by a factor of (only) two. 

A simple prediction of this theory is that that the ratio 
(HJ) does not depend on the sign of the mean magnetic field, 
i.e. it does not vary from cycle to cycle and (for a dipolar 
field) from hemisphere to hemisphere. As a consequence, 
the sign of the cross helicity (u ■ b) should vary from cycle 
to cycle and between the hemispheres. Zhao et al. (2011) 
indeed found indications for a variation from hemisphere 
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to hemisphere in SOHO/MDI magneto-grams and doppler- 
grams recorded in 2000, 2004 and 2007. 



2. Mean-field electrodynamics 

Let U + u and B + b be the fluctuating velocity and mag- 
netic field with the average values U and B. The scalar 
correlation between the fluctuations of flow and field, i.e. 
the cross helicity {u ■ b), is a pseudoscalar. In the same 
sense, the cross correlation tensor (uibj) is a pseudotensor. 
We are here only interested on its symmetric part 



Hij - 



(uibj) + (ujbi) 



As we have shown the tensor can be finite in presence 
of a mean magnetic field B and for density-stratified fluids 
(Riidiger et al. 2011). Consider these quantities as small 
enough that expressions linear in the mean magnetic field 
influence of these quantities are sufficient. The same may 
hold for the shear which influences the (radial) magnetic 
field components. It is then straightforward to formulate 
the relation 

= a(G ■ B)6ij + kCIl, + GJi.i + j(B itj + B jfi )+ 
+a(g-B)(U itj + U jti ) + 
+b{U it iB j + UjjBi)gi + c; / ',./</, + I ,.///, + 
+d(Ui, t Bj + U u Bi)gi + <\f ,.,<,, + I , I h. (3) 

No other formations are possible linear in the mean field B, 
the stratification vector G and the shear of the divergence- 
free mean flow U. For the tensor components we find 



H yz = PgB y + (a + b + c)gB z U y>z 
and 

H zz = (a+ 2/3)gB z + 2egB y U y>z 



(4) 



(5) 



if a box coordinate system (x, y, z) for the latitudinal, az- 
imuthal and vertical direction is introduced. The z-axis is 
aligned with the stratification vector, i.e. represents the ra- 
dial direction in spherical geometry. The x and y coordi- 
nates denote the horizontal directions. Without shear the 
correlation H zz measures the vertical magnetic field while 
the correlation H yz measures the azimuthal field. The cor- 
relations are also influenced by the shear U y . z . With the 
shear included finite values for both the correlations (U]) 
and ([5]) result even for the case that the field has only one 
component. For known values of the correlations, the coef- 
ficients and the vertical field both the azimuthal field and 
the shear can be computed. We cannot, however, be sure 
that all the coefficients a....e must be nonzero. First test 
calculations of H zz under the presencee of horizontal field 
and shear did not yield finite values of e (A. Brandenburg, 
private communication). 

The turbulent flow is assumed anclastic, so that 
div pu — 0. It is convenient to use the Fourier transfor- 
mation of the momentum density m = pu, i.e. 



m(r,t) = J m(k,u) c^ kr ~^dk duj, 

and similarly for the fluctuation of the magnetic field. 



(6) 



The spectral tensor of the momentum density that ac- 
counts for the stratification of the turbulence to the first 
order terms reads 

x (5ij - hkj/k 2 + (Kikj - Kj ki) /(2k 2 )) , (7) 

where k = (z — z')/2, k = z+z', q is the Fourier transform 
of the local spectrum, 



q(k,uj,r) — p 2 E(k,ui,r) — J q(k, u, K,)e 1K r dn, 



(2) so that 



(u 2 ) = J J E(k,uj,r) dkduj. 



(8) 



(9) 



o o 

Derivation of the cross correlation yields 

Hio = \m (G l B J + GjBi) 

3 2 \ 

Y^ T + Y5 ( B 3.* + b <j)i 
where G = Vlogp is the gradient of density and 
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(10) 



(11) 



(12) 



where r\ is the molecular magnetic diffusivity. Both quanti- 
ties remain finite in the high-conductivity limit. 

From the cross correlation tensor (flU|) the cross helicity 
(u-b) = rjT (G ■ B ) is obtained. From Eq. (fTO]) we find the 
slightly more complicated expression 



(u z b z ) = rj T GB z - 



( 3?7t 


h 2f) 






15/ 





(13) 



where G = G z is the only nonzero radial components of the 
density-stratification vectors. Note the negativity of G. An 
upward divergence of the mean field would reduce the effect 
of density stratification but for uniform field components 
the result is (u z b z ) — (u-b). 

A real difference, however, between the both correlation 
expressions is due to a possible gradient G' of the turbu- 
lence intensity it rms . One easily finds that for vertical fields 
the turbulence intensity gradient G' enters the expressions 
for the correlations such as 

(u-b) = (G+±G')mB z , (u z b z ) = ( G +^G') VT B Z .(U) 

In the bulk of the convection zone G' = G' z is positive 
while G is negative. Hence, \(u z b z )\ > \(u ■ b) \ for positive 
B z which is confirmed by the presented simulations (see 
below). 

By elimination of G one finds 



B„. 



(15) 
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The magnetic eddy difFusivity can thus be determined if 
the LHS of (|T5|) is calculated from magnetoconvection sim- 
ulations when the density scale height H p is known from 
numerical models of the solar atmosphere. As only the cor- 
relation (u z b z ) can directly be observed one needs a numer- 
ical model for the application of the LHS of ([T5]) to derive 
the eddy diffusivity at the solar surface. 



3. Numerical simulations 

We perform simulations for a number of different parameter 
combinations. These parameters include the strength of the 
imposed vertical field B z , the viscosity v and the magnetic 
diffusivity coefficient r\. 

The numerical simulations are done using the Nirvana 
code (Ziegler 2002), which uses a conservative finite differ- 
ence scheme in Cartesian coordinates. The code solves the 
equation of motion, the induction equation, and the equa- 
tions of energy and mass conservation. We assume an ideal, 
fully ionized gas that is heated from below and kept at a 
fixed temperature at the top of the simulation box. Periodic 
boundary conditions apply at the horizontal boundaries. A 
homogeneous vertical magnetic field is applied. The upper 
and lower boundaries are impenetrable and stress-free. 

The simulation volume is a rectangular box. The strat- 
ification is along the z- coordinate and it is piecewise poly- 
trophic, with the polytrophic index chosen such that the 
hydrostatic equilibrium state is convectively stable in the 
lower and unstable in the upper half of the simulation box. 
In the following, p denotes gas pressure, p mass density, T 
temperature, g gravity, k thermal conductivity and c p the 
specific heat capacity at constant pressure. 

The gas is initially in hydrostatic equilibrium, i.e. 



Density 



Temperature 



dp 
dz 



pg = o, 



(16) 



where g = const., and the heat flux through the box is 
vertical and constant, 



F 



dT 

-K—— = const. 

dz 



(17) 



The equation of state is that for an ideal gas and the heat 
conductivity is constant in the upper and lower layer, re- 
spectively, but its values differ between the two layers. 

In the dimcnsionless units the size of the simulation box 
is 8 x 8 x 2 in the x, y, and z directions, respectively. The 
numerical resolution is 512 x 512 x 128 grid points. The 
stratification of density, pressure, and temperature is piece- 
wise polytrophic as described in Ziegler (2002). Similar se- 
tups have been used by Cattaneo et al. (1991), Brummell 
et al. (1996), Brandenburg et al. (1996), Chan (2001) and 
Ossendrijver et al. (2001). The initial state is in hydrostatic 
equilibrium but convectively unstable in the upper half of 
the box. The z coordinate is negative in our setup, with 
z = at the upper boundary. The stable layer thus extends 
from z = — 2 to z = — 1, the unstable layer from z = — 1 to 
z = 0. The density varies by a factor 5 over the depth of 
the box, i.e. the density scale height is 1.2. 

Figure Q] shows snapshots of the fluctuations of density 
and temperature for Ra = 10 7 in a horizontal plane close to 
the upper boundary. The density is increased at the bound- 
aries of the convection cells and decreased at the center. 
The opposite is true for the temperature, which is highest 





0.94 0.95 0.96 0.97 0.98 9.846^)1 9.996-01 1.016+00 1.03e+00 1.05e+O0 

Fig. 1. Fluctuations of density and temperature in the 
upper part of the unstable layer at z = —0.05 for Ra=10 7 



at the center of a convection cell and lowest at the bound- 
aries. Vertical velocity is positive, i.e. upwards, at the center 
and negative, i.e. downwards, at the boundaries. The mag- 
netic field is strongly concentrated in a few small patches 
which coincide with cell corners, where the gas horizontal 
flow converges and the vertical flow is downwards. 

The initial magnetic field is vertical and homoge- 
neous. We run the simulations until a quasistationary state 
evolves. Our control parameters are the heat conduction co- 
efficient, k and the Prandtl number, Pr= v / 'k. Convection 
sets in if the Rayleigh number, 



Ra = 



pgcpcfi f dT 



Tkv 



dz 



cp 



(18) 



with the density p, the specific heat capacity cp, the gravity 
force g, and the length scale d, exceeds a critical value. 
The length scale is defined by the depth of the convectively 
unstable layer, i.e. d = 1. After (fTJ) the correlations and 
the mean magnetic field always have opposite signs. This 
has also numerically been realized. For positive values of 
the mean magnetic field B z the cross helicity is negative in 
the unstably stratified layer. If the field polarity is reversed 
and everything else is left unchanged the cross correlation 
becomes positive with the same amplitude. 

The velocity field, which is measured in units of c ac /100, 
shows the asymmetry between upwards and downwards mo- 
tion as it is characteristic of convection in stratified media. 
The downward motion is concentrated at the boundaries 
of the convection cells and particularly at the corners. The 
upwards motion fills the interior of the convection cells (see 
Fig. [2]). As it covers a much larger area the gas motion 
is much slower than in the concentrated downdrafts. The 
magnetic field shows a similar pattern. The vertical field 
is concentrated in the areas with downwards motion and 
weak in the areas with upward motion. As the total ver- 
tical magnetic flux is conserved, this is the result of field 
advection. 

Figures [3] and 0] hold for Ra = 10 7 and for weak and 
strong magnetic fields. The value of both the Prandtl num- 
ber and the magnetic Prandtl number is 0.1. The left dia- 
gram shows the horizontal average of the cross helicity as 
a function of the depth and the right diagram shows the 
same for the correlation of the vertical components, (u z b z ). 
There is a difference between the two quantities, with the 
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Fig. 2. The same as in Fig.Q]but for the fluctuations of the 
vertical flow and the vertical field. 



<u.b> 



<uz.bz> 
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-0.010 



-2.0 -1.5 -1.0 -0.5 0.0 



-2.0 -1.5 -1.0 -0.5 0.0 



Fig. 3. The numerical values for the cross helicity (u ■ b) 
(left) and the coefficient (u z b z ) (right) for weak magnetic 
field, B z = 1CT 3 . 



<u.b> 
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-2.0 -1.5 -1.0 -0.5 0.0 



Fig. 4. The same as in Fig. [3] but for B z = 1. 



vertical component actually being twice the cross helicity. 
Equation (|15|) can thus be written as 



(u z b z ) 
B z 

with 



(19) 



(20) 



The correlations do not vanish abruptly at the bottom of 
the unstable layer because of overshoot, which affects the 



Fig. 5. The same as in Fig.@]but with a reduced numerical 
resolution (128 x 128 x 128 instead of 512 x 512 x 128, left) 
and computed with a later snapshot (right). The differences 
are unsignificant. 

upper half of the stable layer. The correlations are there 
positive and much smaller than in the unstable layer. 

The results in Fig. [3] are given in arbitrary units defined 
by the code. Velocities are given in units of c ac /100 with 
the isothermal speed of sound c ac . With an approximate 
value of c ac ~ 6.6 km/s at the optical depth r = 1 of the 
Sun the simulations lead to the cross correlation velocity 
(u z b z )/B z ~ —10 in units of 0.066 km/s (Fig. [31 right), i.e. 
after (O 



1.1 km/s. 



(21) 



A similar value also holds for magnetic field amplitudes over 
many orders of magnitudes. For much stronger magnetic 
fields Fig. @] (right) yields the slightly smaller value of 0.93 
km/s. 

A characteristic velocity results as the cross correlation 
velocity 



\(vb)\ 

B z 



(22) 



Using the maximal values in Fig. [3J (left) we find U c ~ 
5 in units of c ac /100. Hence, the simulations lead to the 
cross correlation velocity U c ~ 0.33 km/s. The same value 
holds for magnetic field amplitudes over many orders of 
magnitude (Fig. [H left). 

It also makes sense to normalize the cross correlation in 
the form 



(u-b) 



B Z ^W) 



(23) 



which is the ratio of the cross correlation velocity (f2"2"j) and 
the rms velocity of the turbulence. Its numerical value does 
not depend on the internal units of the code so that c v is 
a general and basic result of the simulations. Close to the 
surface the maximal numerical value is ~ 0.5. Test cal- 
culations for various magnetic fields over many orders of 
magnitudes show this value as almost uninfluenced by the 
magnetic-field suppression. Resulting from the overshoot 
phenomenon at the bottom of the unstable layer always 
small negative values there appear. The correlation coeffi- 
cient 



(u-b) 



(24) 
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Fig. 6. The numerical values for u rms (left) and for 6 rms 
(right) for weak magnetic field, B z = (top) and strong 
magnetic field, B z = 1 (bottom). Only the magnetic fluc- 
tuations depend on the background field amplitude. 



for the cross helicity is much smaller than ((23]) as always 



(25) 



very similar to the result of Ossendrijver et al. (2001). The 
relation (|2"4")) proves to be true for all amplitudes of the 
mean magnetic field between 10 -5 and 10. One finds for all 
calculations a characteristic correlation coefficient c ~ 0.1 . 

During the simulations there are significant temporal 
fluctuations. The convective instability initially grows ex- 
ponentially until its saturation when the system settles in a 
statistically steady state. To illustrate the effect of the fluc- 
tuations on the results, Fig.[S]shows {u z b z ) for B z = 1 from 
another snapshot of the same simulation (right) and from a 
simulation with lower resolution. In both cases differences 
for the correlations are visible but rather small compared 
with the numbers in Fig. El 

Figure [6] contains all informations about the kinetic and 
magnetic energies of the magnetoconvection. The rms value 
of the velocity is hardly influenced by the large-scale mag- 
netic field. In physical units we find an averaged value of 
u rms ~ 0.1c ac ~ 0.66 km/s. Contrary, the magnetic energy 
strongly depends on the applied magnetic field. In dimcn- 
sionlcss units it is in both cases b Tms /u rms ~ 0.6B Z which 
leads to 



E 



3600- 



Bl 



kin 



(26) 



in physical units. At the top of the convection zone we find 
very small contributions of the magnetic energy for B z = 1 
Gauss while for 1000 Gauss there is almost equipartition. 



Table 1. Results from the analysis of the SST and 
HINODE data. The resolution is measured on the solar sur- 
face; the SST resolution is reduced by rebinning the data. 



Dataset 


Resolution 


(u z b z ) 


B z 


{u z b z )/B z 






km 


[G km/s] 


[G] 


[km/s] 


[km/s] 


HINODE 


230 


-1.04 


2.55 


-0.41 


0.71 


SST 


115 


-1.82 


2.54 


-0.72 


1.26 


SST* 


172 


-1.57 


2.54 


-0.62 


1.08 


SST* 


258 


-1.32 


2.54 


-0.52 


0.91 


SST* 


343 


-1.05 


2.54 


-0.41 


0.72 



* These data have been rebinned. 



4. Observations 



It is difficult to empirically determine the cross helicity 
(u ■ b) at the solar surface, because it is hard to retrieve 
the horizontal flows and magnetic field components from 
observations. We have, however, the possibility to use the 
relation (u ■ b) « 0.5(u z b z ), known from the above numeri- 
cal simulations. The vertical flow speed and magnetic field 
component can be determined with much better accuracy. 
Then using Eq. (|19|) we can determine the cross helicity 
velocity from the observations. 

For this purpose we have analyses two datasets con- 
taining observations of quiet Sun at disk center, where the 
line-of-sight coincides with the local vertical. Data from 
the CRISP instrument on the Swedish 1-m Solar Telescope 
(SST) cover the 6302.5 A Fel spectral line with 12 equidis- 
tant wavelength positions at 48 mA steps and a contin- 
uum point. They have a pixel scale of 0'.'0592 and a total 
field-of-view of about 60"x60". The second dataset is from 
the spectropolarimcter on the Solar Optical Telescope of 
HINODE and covers both the 6301.5 and 6302.5 A Fel 
lines, has a pixel scale of 0'.'16 and a total (scanned) field- 
of-view of 164" x 328". 

The line-of-sight velocity and magnetic field data for 
the HINODE observations were taken from the level 2 data 
products available onlineQ. Magnetic field strengths have 
been converted to fluxes by taking the filling factor into 
account. The SST data were inverted using the lilia inver- 
sion code (Socas- Navarro 2001). Velocities were calibrated 
using the convective blueshift determined by de La Cruz 
Rodriguez et al. (2011). More details on these two datasets 
can be found in Schnerr & Spruit (2011). 

We show the results for these datasets in Table [TJ The 
cross helicity velocity (u^) as determined from the SST data 
is somewhat higher than that from the HINODE data. At 
least partly this is due to the lower resolution of HINODE 
as compared to the SST. If we rebin the SST data to a 
lower resolution, the cross helicity velocity decreases (see 
Table Q}. The reason for this is that the strongest fields 
and flows arc smoothed out. 



5. Conclusions 

We have shown that nonrotating turbulence at the top of 
the solar convection zone under the influence of a vertical 
magnetic field forms a finite cross helicity. The only con- 
dition is the existence of a vertical stratification of density 
and/or turbulence intensity. The effect would not appear 
within the Boussincsq approximation. It also exists in the 
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high-conductivity limit, i.e. for sufficiently large magnetic 
Reynolds numbers. 

In our understanding the cross helicity is anticorrelated 
to the mean radial magnetic field, i.e. 

(u ■ b) ■ B z < 0. (27) 

For an oscillating dipolar background field the sign of the 
cross helicity differs for both hemispheres and also from 
cycle to cycle. 

The theory can also be used to measure the magnetic 
diffusivity if the cross helicity is known by observations. In 
order to find the cross helicity one has only to correlate 
observed flow fluctuations with observed magnetic fluctua- 
tions. 

The anticorrelation (|2"7| for density-stratified turbu- 
lence has been established by Riidiger et al. (2011) for 
a model of numerically-driven turbulence. In the present 
paper buoyancy-driven magnetoconvection has been simu- 
lated in a box with the NIRVANA code. We find that also 
such a turbulence fulfills the relation (|2"7T) . Independent of 
the applied magnetic field amplitude the correlation coeffi- 
cient ((23)) takes the value of 0.5. The ratio (j25|) of the mag- 
netic fluctuations to the applied magnetic field is always of 
the order five. 

We have also shown that for density-stratified turbu- 
lence the identity (u ■ b) = (u z b z ) holds. So far solar ob- 
servations can only measure the correlation (u z b z ). The 
numerical simulations, however, always lead to the result 
(u z b z ) ~ 2(u- b) so that the observed value of (u z b z ) would 
overestimate the actual cross helicity by a factor of two. The 
reason is the vertical stratification of the turbulence inten- 
sity which at the top of the convection zone is antiparallel 
to the density stratification. Hence, both the correlations 
{u z b z ) and (it • b) are reduced but not by the same amount. 

With (u z b z ) ~ 2(u-b) the value of u v can be computed 
by use of Eq. (|19[). The numerical simulations lead to u n ~ 
1.1 km/s. This result is well confirmed by the observations 
which lead to values between 0.7 km/s (HINODE) and 1.3 
km/s (SST). To estimate the value of the eddy diffusivity 
at the solar surface we shall assume a density scale height 
of 100 km and find values close to 77T — 1-1 x 10 12 cm 2 /s 
for the eddy diffusivity at the surface of the quiet Sun. 
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